function [ phi1,phi2,phi3,phi4 ] = fun_prax3( alfa,beta )

lamda = 19e-2;
k=[cos(alfa)*cos(beta);cos(alfa)*sin(beta);sin(alfa)];

%радиус-вектора АЭ
r1 = [lamda/4; -lamda/4; 0];
r2 = [-lamda/4;-lamda/4;0];
r3 = [-lamda/4;lamda/4;0];
r4 = [lamda/4;lamda/4;0];

%фазовые набеги
phi1 = r1.'*k*(2*pi/lamda);
phi2 = r2.'*k*(2*pi/lamda);
phi3 = r3.'*k*(2*pi/lamda);
phi4 = r4.'*k*(2*pi/lamda);

end

